Dendritic and uniform flux jumps in superconducting films 
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Recent theoretical analysis of spatially-nonuniform modes of the thermomagnetic instability in 
superconductors^ is generalized to the case of a thin film in a perpendicular applied field. We solve 
the thermal diffusion and Maxwell equations taking into account nonlocal electrodynamics in the 
film and its thermal coupling to the substrate. The instability is found to develop in a nonuniform, 
fingering pattern if the background electric field, E, is high and the heat transfer coefficient to 
the substrate, ho, is small. Otherwise, the instability develops in a uniform manner. We find the 
threshold magnetic field, H{iii^{E , ho) , the characteristic finger width, and the instability build-up 
time. Thin films are found to be much more unstable than bulk superconductors, and have a stronger 
tendency for formation of fingering (dendritic) pattern. 

PACS numbers: 74.25.Qt, 74.25.Ha, 68.60.Dv 



I. INTRODUCTION 

The thermomagnetic instability or flux jumping is 
commonly observed at low temperatures in type-II su- 
perconductors with strong pinning. ^"'^ The instability 
arises for two fundamental reasons: (i) motion of mag- 
netic flux releases energy, and hence increases the local 
temperature; (ii) the temperature rise decreases flux pin- 
ning, and hence facilitates the further flux motion. This 
positive feedback can result in thermal runaways and 
global flux redistributions jeopardizing superconducting 
devices. The conventional theory of the thermomagnetic 
instability^'^ considers only "uniform" flux jumps, where 
the flux front is smooth and essentially straight. This 
picture is true for many experimental conditions, how- 
ever, far from all. Numerous magneto-optical studies 
have recently revealed that the thermomagnetic insta- 
bility in superconductors can result in strongly branched 
dendritic flux patterns. ^"^^ 

In a recent paper we examined the problem of flux 
pattern formation in the slab geometry.^ Experimentally, 
however, the dendritic flux patterns are mostly observed 
in thin film superconductors placed in a perpendicular 
magnetic fleld. A flrst analysis of this perpendicular ge- 
ometry was recently published by Aranson et al.'^^ Here 
we present more exact and complete picture of the den- 



dritic instability and analyze the criteria of its realization. 

In the following we restrict ourselves to a conventional 
linear analysis^'^'^^ of the instability and consider the 
space-time development of small perturbations in the 
electric field, E, and temperature, T. In contrast to the 
slab case,^ the heat transfer from the superconductor to a 
substrate as well as the nonlocal electrodynamics in thin 
films are taken into account. Consequently, the results 
depend significantly on the heat transfer rate, /iq, as well 
as on the film thickness, d. Our main result is that the 
instability in the form of narrow fingers perpendicular to 
the background field, E, occurs much easier in thin films 
than in slabs and bulk samples, and the corresponding 
threshold field, Ec, is found to be proportional to the film 
thickness, d. 



II. MODEL AND BASIC EQUATIONS 

Consider the perpendicular geometry shown in Fig. 1, 
with a thin superconducting strip placed in a transverse 
magnetic field, H. The strip is infinite along the y 
axis, and occupies the space from —d/2 to d/2 in the 
z-direction and from to 2u; in the cc-direction. It is as- 
sumed that d-^w. In the unperturbed state the screen- 
ing current flows along the y-axis . The distributions of 
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FIG. 1: A superconductor strip on a substrate (only the left 
half is shown). The dark gray area is the flux-penetrated 
region. 

the current density, j, and magnetic induction, B, in the 
flux penetrated region < x < ^ are determined by the 
Maxwell equation 

curl B = fioi, (1) 

where the common approximation B = /ioH is used. To 
find the electric field and the temperature we use another 
Maxwell equation together with the equation for thermal 
diffusion, 

curlE = -dB/dt, (2) 
C{dT/dt) = kV^T+jE. (3) 

Here C and n are the specific heat and thermal conduc- 
tivity, respectively. 

Equations (l)-(3) should be supplemented by a 
current-voltage relation j — j{E, B,T). For simplicity 
we assume a current- voltage curve of the form 

i=UT)g{E){E/E). (4) 

A strong nonlinearity of the function g{E) leads to for- 
mation of a quasi-static critical state with j « jc{T), 
where jc is the critical current density. ^'^ We neglect any 
B-dependence of jc, i. e., adopt the Bean model. The 
exact form of g{E) is not crucially important, the only 
issue is that it represents a very steep E{j) curve having 
a large logarithmic derivative, 

n{E) = d\nE/d\njKjc/crE:$>l. (5) 

Here a is the differential electrical conductivity, <t{E) = 
dj/dE. The parameter n generalizes the exponent in 
the frequently used power-law relation E oc j" with n 
independent of E. 

The key dimensionless parameter of the model is the 
ratio of thermal and magnetic diffusion coefficients,^ 

T = ^QKa/C . (6) 



The smaller r is, the slower heat diffuses from the per- 
turbation region into the surrounding areas. Hence, one 
can expect that for smaller r: (i) the superconductor 
is more unstable, and (ii) the formation of instability- 
induced nonuniform structures is more favorable. 

In the following we assume that the strip is thinner 
than the London penetration depth. Al, and at the same 
time much wider than the effective penetration length, 

d < Xl -^i Vdw . 

The stationary current and field distributions in a thin 
strip under such conditions were calculated by several 
authors, ^^"^^ finding that the flux penetration depth, i, 
is related to the applied field by the expression 

I / w = TT^H^-/2<fjl . (7) 

Here it is assumed that the penetration is shallow, or 
more precisely that Aeft i -^w. 

III. PERTURBATION ANALYSIS 

A. Linearized dimensionless equations 

We seek solutions of Eqs. (l)-(4) in the form 

T + ST{x,y,z,t), E + SE{x,y,z,t), j + S3{x,y, z,t) 

where T, E and j are background values. The background 
electric field may be created, e. g., by ramping the exter- 
nal magnetic field, and for simplicity we assume it to be 
coordinate independent. Allowing for such a dependence 
would only lead to insignificant numerical corrections, as 
discussed in Ref. 1. Similarly, we will assume a uniform 
background temperature. 

Whereas it follows from symmetry considerations that 
Ex = 0, both components of the perturbation, (5E, will 
in general not vanish. Linearizing the current-voltage 
relation. Eq. (4) one obtains: 

We shall seek perturbations in the form 

ST ~ T* 6 cxp{Xt / to + ikxS, + ikyT]) , 
SEx,y = Ee^^y exp(Ai/to + ikxS, + ikyrf) , 
53x,y = jcix,y exp(At/to + ikx^ + ikyr]) , (9) 

where 9, e and i are ^-dependent dimensionless Fourier 
amplitudes. The coordinates are normalized to the 
adiabatic length a = CT* / ji^ic where T* = 
— {d\Yijc/dT)~^ is the characteristic scale of the tem- 
perature dependence of jc, so that ^ = x/a, rj — y/a, 
( = z/a. The time is normalized to to = <^CT* / j^ = 
fiQcra^, which is the magnetic diffusion time for the length 
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a. Re A is the dimensionlcss instability increment, which 
when positive indicates exponential growth of the per- 
turbation. 

We can now use the formulas (9) to rewrite the basic 
equations in dimensionlcss variables. From Eq. (8) one 
finds for the components of the current density pertur- 
bation i, 

ix = ex, iy = -6 + n~^£y. (10) 

Combining the Maxwell equations (1) and (2), and the 
thermal diffusion equation (3) yields 



k X [k X e] = A?! i , 



\e = T\ -k 



(11) 
(12) 



Magneto-optical imaging shows that flux patterns pro- 
duced by the dendritic instability^"^'' are characterized 
by having ky ^ k^- Therefore, we have neglected the 
heat flow along x direction compared to that along the 
y direction. Later we will check the consistency of this 
assumption by showing that indeed the fastest growing 
perturbation has ky^k^. 



B. Boundary conditions 

We assume that heat exchange between the supercon- 
ducting film and its environment follows the Newton cool- 
ing law. For simplicity we let the boundary condition, 
kV(T + ST) = -ho{T + ST - To), apply to both film 
surfaces. Here Tq and hg are the effective environment 
temperature and heat transfer coefficient, respectively. 
Eqs. (12) and (10) can now be integrated over the film 
thickness to yield 



(1 



n\ + nT{kl + /i) + 1 



where 



2haa^ / nd . 



(13) 



In the remaining part of the paper we let 0, e and i denote 
perturbations averaged over the film thickness. 

We seek a solution of the electrodynamic equations in 
the flux penetrated region, < ^ < l/a. At the film 
edge, ^ = 0, one has Sjx = and, consequently, SE^ = 0. 
In the Meissner state both the electric field and heat 
dissipation are absent, so that SEy = ST = Sjy = at 
the flux front; ^ = £/a. Thus, the Fourier expansions for 
the X and y components of electric field perturbation will 
contain only sin{kxS,) and cos{kx£,), respectively. Then 
the boundary conditions are satisfied for 

kx = {Tra/2£){2s + l), s = 0,1,2,.... 

Since £ depends on magnetic field, the values of k^ are 
also magnetic field dependent. 



Now we can integrate Eq. (11) over the film thickness 
and employ the symmetry of the electrodynamic problem 
with respect to the plane z = 0. It yields 

ikyikx^y ~t- ikySx) —^£x — A?^;^^;, 
-kx{kx£y + ikyEx) + '^£'y = -A7i/(A, fcj;)e;, . (14) 
We have here introduced the function 

^^'"''-Syn n\ + nT{kl + K) + l- 

Note that the equation for the z-componcnt of the field is 
satisfied automatically. The derivatives y with respect 
to C arc taken at the film surface, C = d/2a. To calculate 
them, one needs the electric field distribution outside the 
superconductor, where the flux density is given by the 
Bio-Savart law. 



3 , j X (r - r') 



B(r) = A.oH+£ l^^r' 



The perturbation of flux density is then, 



pt/a poo 

SBx,y = ±MoCrf / de / df^'Gi^ - C, V - ri')S3y,: 

Jo J -oo 



Gi^,il) 



1 



47r ie + + {d/2a)^) 



2^3/2■ 



(15) 



Here we have approximated the average over substitut- 
ing = 0. In this way we omit only terms of the order 
of (d/a)^ ^ 1. The integration over ^' should, in princi- 
ple, cover also the Meissner region, ^' > I /a. Though the 
flux density there remains zero during the development 
of perturbation, the Meissner current will be perturbed 
due to the nonlocal current-field relation. However the 
kernel G(C,'7) decays very fast at distances larger than 
d/a and therefore the Meissner current perturbation pro- 
duces only insignificant numerical corrections. 

The perturbation of magnetic field can be related to 
that of electrical field by Eq. (2), which can be rewritten 
as 



SE'^ y/E ^ ^\nSBy,x/ ^J-oajc 



(16) 



Due to continuity of the magnetic field tangential compo- 
nents Eq. (16) is also valid at the film surface, C, = d/2a. 
Thus it can be substituted into Eqs. (14). The Fourier 
components of the kernel function G{S,,r]) with respect 
to r] can be calculated directly yielding 



GiC.ky) 



kya Ki [ky^e + id/2a)^) 



2ttI 



VeTjd/w 



(17) 



where Ki is the modified Bessel function of the second 
kind. 

The above Fourier expansions in cos(fc2;^) and sm(kx^) 
correspond to the finite interval —2£/a < £, < 2£/a. 
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Therefore we should continue ex,y from < ^ < ^|a to 
this interval and then introduce Go: and Gy as analytical 
continuations oi G{£, — ^' ,ky) having the same symmetry 
as Ex and £y, respectively (see Section VI for details). All 
this allows us to rewrite the set (14) as 

—ikxkyEy + [k^ + XrL)£x 

= {d/2a)XnY,Gx{kx,k'^,ky)ex{k'J , (18) 

K 

{kl + \nf)ey + ikxkySx 

= {d/2a)XnfY,Gy{kx,k',,ky)ey{k'J, (19) 

e/a flja 



Gxikx,K,k )]r 1^ r deGi^-e,ky) 

l^y{lix, Hx: l^y) J Jq Jo 



sin(fc^^) sin(fc;$') 
cos{kx£,) cos{k'^C) 



(20) 



We are interested only in the specific case of very thin 
strip, 



a = d/2l < 1 



(21) 



One can then find analytical expressions for the kernel, 
and it turns out that only its diagonal part, kx ~ k'^, is 
important. 

In this manuscript we present analytical expressions 
up to the first order in a, while the plots are calculated 
up to the second order. The second-order analytical ex- 
pressions can be found in section VI. The kernel (20) can 
be written as 



G X ^y{kx 1 kxj ky^ — ^ 



1 - 7(a, kx)a 



a 



(22) 



where j{a,kx) is a dimensionless function. In what fol- 
lows we shall consider only the main instability mode, 
kx = 7ra/2€, which turns out always to be the most un- 
stable one. For this mode, and in the limit a —^ 0, the 
function "f{a,kx) approaches a constant value « 5. 

Substituting the above expression for G into Eqs. (18) 
and (19) one obtains the dispersion relation for X{kx, ky): 



AiX^ + A2\ + A3 = Q. 



(23) 



Here 



Ai = n-ya , 
^3 = k^T H 



A2 = fcy(l + tAi) + nkl + Ai{hT - 1) , 
x-y +nkl{hT + l/n) + kl{hT-l). 



nkik^.T 



IV. RESULTS 

Let us first consider the simple case of a uniform per- 
turbation, ky = 0. One finds from Eq. (23) that the 
perturbation will grow (Re A > 0) if 
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FIG. 2; Solutions of dispersion equation (23) for small and 
large r, for ct — 0.001 and n=20. 



When the flux penetration region, £, is small, i. e., kx is 
large, the system is stable. As the flux advances, kx de- 
creases, and the system can eventually become unstable. 
The instability will take place, however, only if hr < 1. 
Otherwise the superconducting strip of any width will re- 
main stable no matter how large magnetic field is applied. 
This size-independent stability means that at /ir > 1 the 
heat dissipation due to flux motion is slower than heat 
removal into the substrate. 

Equation (24) further simplifies in the adiabatic limit, 
T ^ 0, when the heat production is much faster than 
heat diffusion within the film or into the substrate. The 
instability then develops at fc^/(7a) < 1, which in dimen- 
sional variables reads as fJ-oj^ld > CT*{tt^ /2^). Assum- 
ing small penetration depth, I <^ w, and using Eq. (7) 
this criterion can be rewritten a,s H > i/adiab, with the 
adiabatic instability field. 



H.. 



'd CT* 



adiab 



W 7/io 



- H: 

w 



;lab 
adiab ' 



(25) 



hr < \ — k^ /7a . 



(24) 



Here -ff^diab ^^'^ adiabatic instability field for the slab 
geometry. '^^ This result coincides up to a numerical 
factor with the adiabatic instability field for a thin strip 
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FIG. 3: Dependences of fcj and fc* on h for n=20, r = 0.01, 
a = 0.001 according to Eq. (26). 



FIG. 4: Stability diagram in the plane electric field - heat 
transfer coefficient according to equation Eq. (28) and condi- 
tion 1 - /ir > for 71 = 30 and a = 0.001. 



found recently in Rcf. 27. 

Solutions of Eq. (23) for perturbations with arbitrary 
ky are presented in Fig. 2. The upper panel shows 
ReA(fcy) curves for r = 0.01 and different values of k^. 
For large kx^ i. e., small magnetic field, Re A is negative 
for all ky. It means that the superconductor is stable. 
However, at small k^, the increment Re A becomes posi- 
tive in some finite range of ky. Hence, some perturbations 
with a spatial structure will start growing. They will have 
the form of fingers of elevated T and E directed perpen- 
dicularly to the flux front. We will call this situation the 
fingering (or dendritic) instability. 

For large r an instability also develops at small kx, 
however in a different manner, see Fig. 2 (lower panel). 
Here the maximal Re A always corresponds to ky = 0. 
Hence, the uniform perturbation will be dominant. The 
uniform growth of perturbations for large r has been re- 
cently predicted in Refs. 1,21 and explained by the pre- 
vailing role of heat diffusion. 

Let us now find the critical k* and k* for the fin- 
gering instability, sec Fig. 2 (upper panel). The k* 
determines the applied magnetic field when the insta- 
bility first takes place, while k* determines its spatial 
scale. These quantities can be found from the require- 
ment max{Re A(fcj,)} = for ky ^ 0. In the limit a <^ 1 
we can put Ai = in Eq. (23) and then rewrite it in the 
form 



A 



-(fc^ + /.)r + ^ 



nk?. 



From this expression we obtain 

fc* = (^V^TTT- /nV^, (26) 
k* = \/nhT + ll\/n + l- y/nhr + 1 j /^ 



The dependences of k* , k* on the heat transfer coefficient 
h arc shown in Fig. 3. One can see that k* is always 



larger than k* implying that fingers of elevated T and E 
are extended in the direction normal to the film edge. For 
/i < 1/r and n > 1 we find k* « n^/'^fc* > fc* « l/^/nr. 
Both k* and k* tend to zero as — > 1/r, while for larger 
h the system is always stable due to fast heat removal to 
the substrate. It follows from Fig. 2 that for large enough 
r the instability will develop uniformly, while for small 
T it will acquire a spatially-nonuniform structure. Let 
us find now the critical value Tc that separates these two 
regimes. It can be obtained from the equality Re X{kx = 
k*,ky ^ 0) = 0. When it is fulfilled ReA = both for 
ky = and for ky = k* 0. We find using Eq. (23) that 
the instability will evolve in a spatially-nonuniform way 
if 



T < Tc 



1 



fc*V7«) /h. 



(27) 



Substituting here a and k* we find a transcendental re- 
lation between Tc and h. For n ^ 1 it reduces to 



^1 + VTlTc^ = TTa/^d . 



(28) 



Using this result we can construct a stability diagram in 
the E — ho plane shown in Fig. 4. The curved line marks 
the critical electric field Ec{ho) that separates two types 
of instability: fingering {E > Ec) and uniform [E < Ec). 
This line is calculated from Eq. (28), where the electric 
field is expressed via t a,s E = jcHon/nCT according 
to Eqs. (5) and (6). The straight line is given by the 
condition hr ~ 1. Below this line the superconductor is 
always stable, as follows from Eq. (24) for the uniform 
perturbations, and from Eq. (26) for the nonuniform case. 
At a certain value ~ /icrit, the two lines intercept. We 
find 



^crit — 



and the critical electric field Ec for ho 



is 



Ec{0) 



(29) 



(30) 
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FIG. 5: Stability diagram in the H-E plane according to 
Eq. (24) and Eq. (26). 



while E,{h,nt) ^iEciO). 

For any point {ho,E) belonging to the stable phase 
in the stability diagram, Fig. 4, the flux distribution is 
stable for any applied magnetic field. For the points 
belonging to unstable phases, the instability develops 
above some threshold magnetic field, either Hfing{ho, E) 
or i7uni(/io, E) for fingering or uniform instability, respec- 
tively. Shown in Fig. 5 are three sets of Hfing{E) and 
Huni{E) curves for different values of h^. They represent 
boundaries between the three phases, stable and unsta- 
ble with respect to either fingering or uniform instability, 
as shown in the inset. Using Eq. (7) one can rewrite the 
expression Eq. (24) for i/uni as 



adiab 



1 - 



2T*h 



jidjcE 



-1/2 



(31) 



In the absence of heat removal to the substrate, /lo = 0, 
we obtain the adiabatic instability field, Eq. (25), and 
the iJuni(-E') curve becomes a horizontal line.^* 

The threshold magnetic field for the fingering insta- 
bility, -fffing, is calculated from Eq. (26). A simplified 
expression obtained for /i <C 1/t and n 3> 1, 



jcd j K T*jc 



1/2 



(32) 



shows that at large electric fields iJfing decays as E~^^'^. 
At ho < /icrit the curves Hfij^g{E) and iJuni(-E') intercept 
at the critical electric field Ec determined by Eq. (27). 
At ho > /icrit we have Hfing{E) < Huni{E) for any E, so 
the lines do not intercept and the instability will develop 
into a fingering pattern. 



V. DISCUSSION 

Let us compare the present results for a thin film in 
a perpendicular magnetic field with results of Ref. 1 for 



a bulk superconductor. In both cases the instability de- 
velops into a fingering pattern if the background electric 
field in the superconductor exceeds some critical value 
Ec- The values of Ec are however different. Their ratio 
for a thin strip and a slab. 



2 



E: 



slab 



CT* 



(33) 



is expected to be much less than unity. For jc = 10^*^ 
A/m2, C = 10^ J/Km3, n = IQ-^ W/Km, T* = 10 
K and d = 0.3 /im, we find from Eq. (30) that Ec ~ 
4 • 10""* V/m, while according to Ref. 1, E^^^^ = 0.1 
V/m. Consequently, the development of thermomagnetic 
instability into a fingering pattern is much more probable 
in thin films than in bulk superconductors. 

The threshold magnetic field for the fingering instabil- 
ity, fffing, is also much smaller for thin films. Comparing 
Eq. (32) with the results of Ref. 1 for a slab^ we find 



fing 



V2 d 



;lab 



(34) 



Here I* = {TT/2)y^KT* /jcE is the flux penetration depth 
at the threshold of the fingering instability, H = Effing- 
Experimentally, the fingering instability always starts af- 
ter the flux has penetrated a noticeable distance from 
the edges, such that I* 3> d.^"^" Hence, for a thin film 
the fingering instability should start at much smaller ap- 
plied fields than in bulk samples (by a factor of ~ 10^ 
for films with d ~ 10~*w). The difference between 
the threshold fields for the two geometries here is even 
stronger than for the case of uniform instability in the 
adiabatic limit, see Eq. (25). Assuming the above val- 
ues of parameters and w = 2 mm we find from Eq. (32) 
that Hfing[Ec{0)] = iJadiab ~ 1 niT. This value becomes 
larger if we take into account the heat transfer to the 
substrate. It is therefore in excellent agreement with 
experiment, ^"^^'^^'-'^^'^^'^^'^^ where the threshold field is 
typically of the order of a few milliTesla. 

The spatial structure of the instability predicted by our 
linear analysis is a periodic array of fingers perpendicu- 
lar to the film edge. Its period can be estimated from 
Eq. (26). For E ^ Ec, h ^ and n > 1 one finds 



TT^CT* 



2jn^/^^ioPcd 



(35) 



which yields dy w 100 /im for n ~ 30. Numerical analysis 
of the instability development shows^'^^ that beyond the 
linear regime the periodic structure is destroyed and only 
one (strongest) finger invades the Meissner region. This 
scenario is indeed reproduced experimentally, and the ob- 
served width of individual fingers, 20-50 /im,^'^'*^'^^ is 
very close to our estimate, dy/2. 

The finger width and the threshold magnetic field also 
depend on the dimensionless parameter h characteriz- 
ing the thermal coupling to the substrate, Eq. (13). In 
turn, h, grows rapidly with temperature because of a 
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strong T dependence of C and jc- One can therefore 
make several testable predictions from the dependences 
and k*(h) shown in Fig. 3: (i) There must be a 
threshold temperature Tth above which the instability is 
not observed, (ii) When approaching Tth, the instabil- 
ity field diverges since fc* — > 0. (iii) When approaching 
Tth the characteristic width of individual fingers increases 
since k* — > 0. The last prediction has also been obtained 
in the boundary layer model allowing calculation of the 
exact finger shape. The first and the second predic- 
tions have already been confirmed experimentally. ^"'^^ 
As for the last one, the T dependence of the finger width 
has not yet been studied. At the same time, there is 
a solid experimental evidence®'^°'^^'^° for an enhanced 
degree of branching as T — > Tth that can be quantita- 
tively described as a larger fractal dimension of the fiux 
pattern.^" This abundant branching could be an indirect 
consequence of the increased finger width since wider fin- 
gers are presumably more likely to undergo splitting. 

The present problem of fingering instability in a thin 
film has two new features compared to a similar problem 
for a bulk superconductor, (i) nonlocal electrodynamics 
and (ii) thermal coupling to the substrate. The non- 
locality results in much smaller values of the threshold 
magnetic field i/flng and the critical electric field Ec in 
films than in bulks. If a film is made thinner, it becomes 
even more unstable since iJfing oc d, and has a stronger 
tendency to form a fingering pattern since £'c(0) oc dP. 
The thermal coupling to the substrate has a somewhat 
opposite effect. It can lead to an ultimate stability if 
h > l/r - a situation that is never realized in bulks. A 
moderate coupling, /i <C l/r, slightly renormalizes -ffgng 
and Ec, i.e. makes the film a little bit more stable and 
less inclined to fingering. 

Let us now compare the results presented in this work 
to those obtained in a similar model by Aranson et 
al.'^^ Our expressions for the "fingering" threshold field , 
Eq. (32), and for the finger width, Eq. (35), agree with 
their results up to a numerical factor. For r 3> 1 our 
results for the "uniform" threshold field (derived from 
Eq. (24)) are also similar to results of Ref. 21. As a 
new result, we find that there exists a critical value of 
the parameter r, Eq. (27), which controls whether the 
instability evolves either in the uniform, or in the finger- 
ing way. Shown in Fig. 4 is the stability diagram where 
the line Ec{ho) separates regimes of fingering and uni- 
form instability. Other new results of this paper are: (i) 



the existence of a field- independent "critical point" , ft-crit, 
such that for /iq > /icrit the instability always develops 
into a fingering pattern, and (ii) the full stability diagram 
in the H-E plane. Fig. 5, containing all three phases. 

The background electric field needed to nucleate the 
fingering instability can be induced by ramping the mag- 
netic field, E ^ HI (X HH^ for I ^ w, where H is the 
ramp rate. This is the lowest estimate since the fiux pen- 
etration in practice is strongly nonuniform in space and 
in time,^^ and there can be additional sources of E due 
to random fluctuations of superconducting parameters. 
The occurrence of the fingering instability even at rather 
low ramp rates^'^"^^'^^'^^"^" is therefore not surprising. 

The build-up time of the instability can be estimated as 
to ~ 0.1/zs if the flux-flow conductivity a = 10^ il~^m~^. 
Our linear analysis assumes that the perturbations of T 
and E grow in amplitude, but remain localized within the 
initial flux penetrated region. Numerical results show^'^^ 
that at t ^ to the perturbations also propagate into the 
Meissner region. This propagation can be described by 
recent models^^''^^ that predict a characteristic propaga- 
tion speed in agreement with experimental values of 10- 
100 km/aJ'^ 

111 conclusion, the linear analysis of thermal diffusion 
and Maxwell equations shows that a thermomagnetic in- 
stability in a superconducting film may result in either 
uniform or finger-like distributions oiT, E and B. The 
fingering distributions will be observed if the background 
electric field E > Ec, where Ec grows with the film thick- 
ness, the critical current density, the thermal conductiv- 
ity and the thermal coupling to the substrate. Due to 
nonlocal electrodynamics in thin films they turn out to 
be more unstable than bulk superconductors and more 
susceptible to formation of a fingering pattern. 
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VI. APPENDIX 

Here we present in more detail some derivations omitted in main text. 



A. Eq. (11) 



In components, equation (11) reads as 



iky ~ ikye^j - = -Xue^, , (36) 



d f de 



g^yg^ ■ ^kySy] = 0. (38) 



B. Transition from Eq.(17) to Eq.(20) 



The Fourier transformation for the kernel function G with respect to variable rj reads 



G{tv,a) = / ^G{^,ky, a)e^'y\ (39) 



Gi^,ky,a) = / d7^G{^,v,a)e-'''y\ (40) 



Substituting Fourier transformation of kernel G into Eqs. (15) and (16), we find 



e;(e, ky) = 2a'Xn f''' d£,'G{i - C, ky)e,iC, ky), (41) 





4(e, ky) = 2a^\nf{\ ky) [ <G(^ - ky)ey{^\ ky). (42) 

Jo 

These are two linear independent integral equations with difference kernel. Each of these equations can be solved 
separately by a standard method. 

Although the procedure of solution of Eqs. (41) and (42) is standard, we describe it in details to avoid any 
mistakes. In a strict sense, the symmetry of the Fourier expansions with cos(fc2;^) and sin(fc2:^) corresponds to the 
interval —2i/a<£^<2£/a for these trigonometric functions to be orthogonal at different k^- Thus, we should formally 
continue analytically Sx^y on this interval. These continuations are different for the sine and cosine functions. It could 
be easily found that the symmetry of the function ej,(^) corresponds to the following continuation ej,(^) from the 
interval < ^ < ^/a to ~2i/a < ^ < 2£/a 

£^(^) = ey{2e/a ~ at £/a<C< 2e/a (43) 

and ej,(^) is an odd function with respect to ^ = 0. The continuation s^xiO of the function ex(C) is 

£x(0 = ~exi2e/a - ^ at £ < ^ < 21 /a (44) 

and is even with respect to ^ = 0. The functions Sx.yi^,) coincide with Sx,y{^) at < ^ < £/a and can be expanded 
using the same trigonometrical functions. So, we get 



1 



2 



2i/a 



ex{kx) = T; ex{0cosikx0d^^2 e,(0 cos(fc,e)rfe, (45) 



2£/a Jo 
■2e/a c-lja 



^vikx) ^ t; £yiO ^^kxOd^ = 2 ey{^)sm{kxi)d^. 

^ J -2l/a Jo 



(46) 
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Let us produce a similar continuations with the function G. First, we consider G{£^ — in Eq. (41) as a function 
of two independent arguments, G(^,^'). Then we can define the continuation Gx{C) of as 

G^iC) = -G{2£/a - at £/a < < 2£/a (47) 

and even with respect to = 0. Second, in Eq. (42), we can define the continuation Gy{(,') of G(^') as 

Gy{^') = G{2e/a - C') at l/a < ^' < 2£/a (48) 

and Gy{^') is an odd function with respect to ^' = 0. The functions Gx.y{S.,C) obey the same symmetry with respect 
to the other argument ^ since G[(^ — £,')] is symmetric with respect to a permutation of ^ and ^' as follows from (15). 
The integrals, Eqs. (41) and (42), can now be written as 

a'^Xn f^^'" 

e.'iO^^ dC'G.(C,m(e'), (49) 

sy'io = ""'^"^'^^^ deGy{^,aeyia m 

^ J-2l/a 

The solutions of these equations coincide with the solutions of the original Eqs. (41) and (42) at < ^ < ^/a and 
have the symmetry necessary for the suggested Fourier transformations. 
Let us introduce the Fourier transformations 



G.iLO = E G,(fc„fci)cos(fc,C)cos(fci^'), (51) 
GyiLO = E G,(fc„fc;)sin(fc.Osin(fcU'), (52) 



where 



G,(fc„ fc;) = 4 f'^ f'^ dCG,{^, a cos(fc,0 cos(fcie'), (53) 
Jo Jo 

Gy{kx, K)=^[' [' d^'Gyi^, O sin(fc,e) MKO, (54) 
Jo Jo 

and Ga:^y{kx,K) = Gx,y{K,kx)- Equations (53), (54) lead to (20) 



C. Properties of kernel G 

First, note that we do not disregard the value of the order of [d/tY i"^ the denominator of the function G. It seems 
as an exceeding of accuracy since previously we neglect the z dependence of electric field and temperature across the 
film. In other words, in the film 

£«(C)=£.(0) + ^^ + ..., (55) 

where the first derivative is omitted due to the symmetry. We neglect all the terms in the right hand side of Eq. (55), 
except the first. That is, we should seemingly neglect all the values of the order of ((i/£)^ and for the temperature 
perturbation as well. However, it is not so since the derivatives across the film thickness includes additional smallness. 
In the case of electromagnetic values it is due to the assumption that d < and for the temperature due to the 
smallness of Bio number Bi = h^d/n in any realistic situation for a thin film. 

As in the previous study^, we analyze here the stability in the linear approximation following two different ways. 
In qualitative approach, we do not specify the boundary conditions exactly assuming that kx is of the order of a/i, 
and analyze evolution of perturbations with each k^ independently. In this approach we evidently need only diagonal 
component of the kernel G with kx = k'^. If we specify the boundary conditions and try to find exact expressions for 
the stability criterion and characteristic time and spatial scales, we need to know kernel function in a more general 
form. Both the approaches should give rise to qualitatively the same results. 
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D. Transition from Eqs. (18) and (20) to formula (22) 

Performing the integration of Fourier transformation for the kernel function G we get 

G{L ky) = (hVeT^) , (56) 

where Ki is the modified Bessel function of the second kind. In general case, the kernel function G{kx, k'^, ky) could 
not be found in the explicit form and Eqs. (18) and (19) should be solved numerically. However, we are interested 
here only in the specific case a ^ 1. Within these limit the analytical expressions for the kernel can be found. 
The expressions for the Fourier components of the function G read 



2fc.„a f'/'' f"'' ,^1 ky,/W^W+^ cos(fc,Ocos(fc;e') 
G,{K,K,ky) = ^ dU de ^ J, , , (57) 



Jo Jo + ' 



Gyik.X,ky)^^J^ d,J^ ^^=^ . (58) 

Under conditions specified for (57), (58), the main contribution to these integrals evidently comes from the region 
1^ — ^'1 < a. At larger |^ — the functions under integrals decays exponentially. So, we can replace the Bessel 
function Ki{x) under integrals by its expression at small value of argument x up to the order of x 

Ki(x) - + ^\nx (59) 
X 2 

From Eqs. (57) and (58) we find in the main approximation accounting for the first term in the rhs of Eq. (59) 

- /•'^°'%. /os(fc,g)cosK (e + t^)] 

2a 1'^'" ^^j""'^ ,..Mk.O Mk'Ai + u)] 



Giik.X,ky) = d^l (60) 



where we proceed to integration over u = ^' — S,- 

We need to compute following integrals, the first one for G!J, and the second for Gyi 

fi/o. re/a-i COS k'u /■^/° /-^/i-C sinfc'w 

/ cos kxS, cos k'^^d^ / du— — / cos fcj,^ sin fc^^d^ / du— — ^ (62) 
Jo J^ u + a Jq u + a 

/ sinkx£,smk'£,d^ / du— — ^+ / sm k^S, cos k fd£, / du— — ^ (63) 
Jo Je " + " "'o Je " + " 

The calculation of this integrals are almost identical, so we show only computing of the first one. In the begining we 
need to compute following integral 

cosfc^rw f'^^"-^ du [''"^^ , 1-cosfc^u 

du -^—^ = / 7^—? - / du ———— . (64) 



^ + J _^ + oi^ J _^ + a- 

The first integral yields [arctan(^/a) + arctan(£/a ~ S_)/a]. Thus we have to calculate 

2a /■^/'' 

— / cosfcjj^cosfc^^ [arctan(^/a) + arctan(i'/a — ^)/a] . 

"■^ Jo 

One can show that off-diagonal integrals, i. e. for kx ^ k'^, are very small, and we will keep only diagonal elements 
which are the same for all kx- Let us denote them as [1 — c(a)]/a. For a ^ 1 one finds c(a) oc a, while a good 
approximation for a ^ 1 case is c(a) w 1.6a°'®*. Even though c(a) — > for a — !■ its account is important since c{a) 
will enter the dispersion law with the factor oc fc^. 
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FIG. 6: Plot of the product 7a versus a for the lowest instability mode, kx — na/2£. Shown in inset is the plot for small values 
of a. 



The second integral in Eq. (64) can be calculated numerically, it has the order of k^, let us denote this integral as 
r(a, kx). The values of r(a, k^) in the limit a <C 1 depend on k^, namely r(a, k^) = for fc^; = and r(a, k^) — > k^ 
for kx ^ 1- 

To complete computation of Eqs. (60) and (61) we also need to calculate integral 

cos fcj;^ sin fc^^d^ / du— ^ — (65) 



which for k.^. = (7ra/2£)(2n + 1) is equal to zero. The same results one can obtain for G°. Despite the fact that Eq. 
(60) contains cos and Eq. (61) contains sin the final results are the same thus we employ the approximation 



^x,yi^xj k,jr. \ ky) 



1 - c{a) 



r{a, k^) 



(66) 



Denoting c{a)/a + r(a, k^) as 7(0:, kx) we get equation (22) in the main text. The dependence of 7(a, k^) on a for 
the main instability mode kx = 7ra/2£ is shown on Figure 6. 



E. Calculation of G up to the second order 

To calculate we need to take following integral: 



/ cos kx£,d£, 



dC Hky VW^W+^) COS kW 



(67) 



To do this we have to spht ln(fc.y\/(^ — ^')2 + a^) into In(fcj^) + ^ ln((^ — ^')^ + a^). The first part of integral can 
be solved analytically quite easily and using that kx = (1 + 2n) we get following result: 



kyU In ky 
■Kikl 



(68) 



The second part of integral (67) cannot be solved analytically. Using numerical computation for a ^ 1 we find an 
approximately linear dependence of this integral on a: a + ha. For kx = (7ra)(2£) this will be —1.386 + 3.1a. The full 
equation for Gx,y will be: 



^a:,2/(^2:; kx^ ky) 



a Trfcf 27r 



(69) 



13 



Using (69) we can get more precise equations for (23). The changes wih be very simple - ah terms containing 70; will 

In k 

transform into (70; — a "^^2 " — '^2Tr^'^ ba)): 



kn, In kti k„, 

A2 = kl{l+TAi)+nkl + Ai{hT-l), 

Az = kyT + nklklT + nkl{hT + l/n) + kl{hT 



F. Exact expression for hcrit 



''crit 



7r2T*2fc2C2 



[n + y/v? + n + 0.5 + V7 • 10"^"n) (71) 



